pub fn new_birthday_probability(n: u32) -> f64 {
    const N: f64 = 365.;
    // p = 1 - N!/((N-n)!N^n)
    1. - (0..n).map(|i| (N - i as f64) / N).product::<f64>()
}
